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Abstract 

Oblique propagation and head-on collisions of solitary structures is studied in a dense mag- 
netized plasma comprised of relativistic ultra-cold electrons and positrons and positive dynamic 
ions using conventional extended multi-scales technique, in the ground of quantum hydrodynamics 
model. The variations of head-on collision phase-shift as well as the characteristic soliton ampli- 
tude and width is evaluated numerically in terms of other plasma parameters such as mass-density, 
normalized magnetic field strength, its angle with respect to the soliton propagation and the rel- 
ative positron number-density. The relevance of current investigations, with appropriate plasma 
parameters for the astrophysical dense magnetized objects such as white-dwarfs, is addressed. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. BACKGROUND 



Extreme conditions in plasma such as high temperatures and pressures can lead to dif- 
ferent features of nonlinear wave dynamics, including propagation and collision, from those 
encountered in ordinary ionized environments 1] . It is well known that dispersive dynamics 
arise due to quantum diffraction phenomenon in many of degenerate plasma kinds [2- 4] . 
It has also been confirmed that in a degenerate plasma there may be critical values which 
govern the nonlinear dynamics [5j. The degeneracy feature which is a fundamental aspect 
of ordinary solids arise due to exclusion mechanism when the de Broglie thermal wave- 
length \b = h/ {2imi e kBT) l l 2 is comparable or higher than inter-particle distances [6[ and 
leads to a much higher pressure on degenerated species compared to the non-degenerated 
ones. It has been shown j3] that in a completely degenerate astrophysical object such as 
a white dwarf the electrons, which follow Fermi-Dirac statistics rather than that of Boltz- 
mann, become relativistically degenerate due to gigantic inward gravitational pressure and 
this can lead to softening of the degeneracy pressure giving rise to the ultimate collapse of 
the star. The effects of the relativistic degeneracy has recently been considered in electron- 
positron-ion plasmas on large-amplitude nonlinear wave dynamics by number of authors in 
ultra- relativistic and non-relativistic degeneracy limits js, 9]. 

On the other hand, more recently, the dynamics of nonlinear ion waves has been investi- 
gated in electron-ion plasma considering a wide range of relativistic degeneracy and it has 
been confirmed that the relativity parameter plays a crucial role in propagation of such 
waves in dense Fermi-Dirac plasmas 10J. However, to our current knowledge, effect 
of the relativistic degeneracy parameter, R (defined in next section), which is 
directly related to quantities such as the electron number-density, n, and mass- 
density of plasma, p, has not been considered in collision dynamics of solitary 
structure. Historically, discovery of the remarkable shape-preservation feature of these 
waves during interactions by Zabusky et.al. [ll| in 1965, made their first important appli- 
cations in communication technology. Later, Washimi et.al 12] showed that such waves, in 
a weakly nonlinear regime, can be mathematically modeled by the well known Korteweg-de 

n 

Vries (KdV) equations and Oikawa et.al. [13[ have used an extended approach to consider 
the interaction of such waves as the superposition of two single KdV-type solitons. 

Electron-positron-ion plasmas are important because they appear in many astrophysical 
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environment such as, active-galactic-nuclei, pulsar magnetospheres, neutron stars and super- 



novas etc. 



14H18j . They may even play important role in evolution of the early Universe 



19] . Electron-positron- ion plasmas can also be produced in intense-laser matter interaction 



experiments 
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J-|23|. It has been noted [l| that the characteristic nonlinear wave frequencies 
are comparably higher compared to that of pair annihilation in such dense plasmas so that 
these waves survive in these extreme environments. In the current study we consider 
in a magnetized electron-positron-ion plasma the propagation and head-on collisions of ion- 
acoustic solitary waves in a wide range of relativity parameter, i.e. plasma mass-density, 
thus the findings may well be applicable to inertial confinement fusion electron-positron-ion 
plasmas. The article is organized in the following way. The basic normalized hydrodynamics 
equations are introduced in section [TT1 Evolution equations along with collision parameters 
are derived in section IIHI Numerical analysis and discussion are presented in section IIVI 
and final remarks are drawn in section |V] 



II. QUANTUM HYDRODYNAMIC FLUID MODEL 

Considering a quantum fully degenerate plasma with relativistically degenerate elec- 
trons/positrons, we will use the conventional quantum hydrodynamic (QHD) fluid equations 
to describe the dynamics of nonlinear excitations. The pair-annihilation rate is ignored com- 



pared to characteristic plasma frequencies [l] and the plasma is considered as collision-less 
due to Pauli-blocking mechanism. Therefore, the closed set of QHD equations may be 
written in the following dimensional form 

fa + V • (fiiVt) = 0, Vi = mi + jvi + kwi, 
^ + V • (rijVj) = 0, Vj = iuj +jvj + kwj, 

^ + (v, ■ v) v, = -^w - gf fe) 7 ~ 2 ^ + 5 o(V, x k), (1) 
- (^r + ( v r V) v/) = -^V0 - -J-VP, + ^(&) v(^p-) , 



V 2 = -47re s i n j + 



where, B = -E>ok, j = {e,p} and Sj = {— 1,+1}, for electrons and positrons, respectively, 
and h is the normalized Plank constant. It is noted that, in fully-degenerate configuration 
from Fermi-Dirac statistical model, the relativistic degeneracy pressure is expressed in the 
following general form, which is valid to the electrons/positrons with arbitrary degree of 



relativistic degeneracy [24 ] 
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TTTTl-C ( l r I [ 1 

Pj = ~dr \ Rj ^ - 3 ) v 1 + ^ + 3 ln n + v 1 + R2 A ) ' (2) 

in which the relativity parameter Rj = [nj/no) 1 / 3 [? ] (n = 87r ^g c ~ 5.9 x 10 29 cm~ 3 ) 
is the ratio of electron/positron Fermi relativistic momentum ppj to m e c. Thus, 
with the new relativistic degeneracy pressure term the basic equations may be rewritten as 

drij 
dt 

^ + v • K-v,.) = o, 

,7-2 



ft + V • (fiiVO = 0, 



^ + (V, • V) V, = -^V0 - (^) (^) 7 " Vn, + 5 (V 4 x k), 
S + CVi • V) V,) = -^V0 - V [l + Rlafn^] ^ + 2} (f ) V 




(3) 



j 

where, i?o = (^eo/^o) 1 ^ 3 (otj = {l,a},a = n p o/n e0 for electrons and positrons, re- 
spectively) is a measure of the relativistic effects (called the relativity parameter) and 
= h^Upi/lrrieC 2 ) 2 , is the quantum diffraction coefficient which is related to the rel- 
ativity parameter via (m e /mi)H^/2 ~ 4.6 x 10 _10 i?o. On the other hand, the parameter 
Rq is related also to the mass-density (of white dwarf, for instance) through the relation 
p ~ 2m p n e0 (l — a) or p(gr/cm 3 ) = (1 — a)p -R(j with p (gr / cm 3 ) ~ 1.97 x 10 6 , where, m p 
is the proton mass. Note that, the density po is exactly within the range of mass-density 
of a typical white dwarf and a contribution of (1 — a) is included in mass-density definition 
for the electron-positron pair production/anihilation effect. The density of typical white 
dwarfs can be in the range 10 5 < p(gr/cm 3 ) < 10 8 , which, neglecting electron-positron pair 
production/anihilation, results in values of 0.37 < Rq < 8 for the relativity parameter. 

Before proceeding with calculations, a clear definition of the relativistic de- 
generacy and distinction between a low-pressure relativistic plasma from the rel- 
ativistically degenerate quantum Fermi-gas is in order. The relativistic degener- 
acy is a completely quantum phenomenon ruled by the uncertainty principle and 
is increased due to the decrease in inter-fermion distances in degenerated plas- 
mas. Although the relativistic effects arise due to increase in fermion number- 
density in a dense degenerate plasma, however, unlike for the low-pressure rel- 
ativistic plasmas the degeneracy pressure in the fermion momentum fluid equa- 
tion usually dominates the relativistic dynamic effects in super-dense plasma 



state. Chandrasekhar [7], combining the relativity and the quantum statis- 
tics, showed that for dense degenerate Fermi-gas such as a white-dwarf with a 
mass-density, p, the degeneracy pressure turns from Pj oc p 5 / 3 (with polytropic 
index 3) dependence for normal degeneracy for the limit R — > to P d oc p 4//3 
(with polytropic index 2/3) dependence for ultra-relativistic degeneracy case in 
the limit Rq — > oo. The relativistic degeneracy starts at mass density of about 
4.19 x 10 6 (gr/cm 3 ) of the order in the core a O.3M white dwarf, which corre- 
sponds to a Fermi- momentum Pp e ~ 1.29m e c (corresponding to the relativistic 
degeneracy parameter value of R ~ 1.29) or the threshold velocity of up e ~ 0.63c 
(the Fermi relativistic factor ^ Fe ~ 1.287). Now comparing the terms in the 
momentum equation Eq. fl3]) reveals that the term containing mass ratio is 
still negligible compared to the degeneracy pressure term. The similar treat- 
ment of ultrarelativistic quantum plasma as in our case has been considered 



in Refs. 



y 



251 ]. However, this assumption breaks down for a relativistic low- 



pressure plasma and one has to follow a coherent relativistic formalism given 
in Refs. [26l. [27]] . Thus, for the range of the relativistic degeneracy parameter 
-Rq ~ p/po < 10, which is used in our analysis (see discussion), it is safe to neglect 
the relativistic dynamic effects compared to that of the dominant relativistic 
degeneracy pressure. On the other hand, it should be noted that the current 
analysis includes the whole range of degeneracy parameter, R , and is not limited 
only to the relativistic degeneracy range. 

The normalized set of QHD equations can be obtained from the following scaling 



~ 2 _ 

V —V, t — , n n n, V c sr V, <p '-^—(p. (4) 



t _ m e c 



where, c, oo P i = \jAite 2 n e o/mi and c sr = \J m e c 2 /mi are the vacuum light speed and the 
characteristic plasma frequency and relativistic sound-speed, respectively, and the parameter 
n e o denotes the electron equilibrium density. Now, neglecting the terms containing the small 
mass-ratio, m e /mj and assuming that fc^Tj m e c 2 , we arrive at the following simplified set 
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of basic normalized equations 



drij 

at 



+ v • (mVi) = o, 

^ + V • (n/V,) = 0, 
"v, + ( V . . V ) V, = -V0 + cu(V, x k) 



at 



0. 



1 + 



1/2 



(5) 



71,; 



where, ZJ = uj C ijuj vi with the ion-cyclotron frequency of u C i = eB /rrii. Solving Eqs. (jSJ) 
for the electron and positron number-densities in terms of electrostatic potential with the 
appropriate boundary requirements ( lim rij = a* and lim = 0), results in the following 
energy relations 



n 



R 3 



3/2 



n p = a R 



3/2 



R 2 ^ 3 _ ^ 2y /i + R 2 a 4 / 3 _ 

where, the equilibrium charge neutrality condition is given by Poisson's relation as 
a + (3 = 1, a = — — , /3 



(6) 



n i0 
n e o' 



(7) 



Transforming the normalized plasma equations (Eqs. (J3])) to the appropriate strained coordi- 
nate defined below admits successful separation of variables and allows elimination of secular 



terms, which consequent 



collision parameters 
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v leads to the desired evolution equations and the corresponding 
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£ = e(kx + ly + mz - c^t) + e 2 P (r], r) + £ 3 Pi(£, r], r) + . . . , 
r\ = e(kx + ly + mz - c v t) + e 2 Q (£, r) + e 3 Q 1 (£, rj, r) + . . . , 



(8) 



T 



e 3 t, c ? 



C, Cn 



where, the interrelated functions Pi and Qi (I = 0, 1, 2, ...) denote the phase of the traveling 
solitary waves to be determined together with the evolution equations in next section. The 
solitons are assumed to travel at directions described by cosine indices (k,l,m) and the 
direction of magnetic field with respect to the collision line is presented by 7 angle which is 
determined in the following relation 



m = cos 7, 

k 2 + I 2 + m 2 = 1. 



(9) 
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The reduced set of plasma equations in new strained coordinate is presented in appendix A. 
From the lowest-orders in e we obtain the following relations 



c| -| + l?)"' 1, + m(1 - a) (l + 4 U ' !1, = ' 



(lla) 



nf } = 3i? " 2 



l + ^ + a -l/3^ 1 + jR 2 a 4/3 

which leads to the following first-order approximations of plasma variables 



(i) 



3R 2 



1 + R2 + a-V3^1 + jR 2 a 4/3j [ ¥ ,(l) (e> r) + ^(l) (r?; r) ] j 

^ = "I [^W(^r) + ^ 1 )(r ? ,r)] , = | [%>«(^ r) + %>«(77, r)] 
(i) _ 



:3i?, 



-2 



/ TTi?I + a- 1/3 v /l + ^o« 4/3 [<^ (1) (^)-y? (1) (^ 



m(l-a) 

On the other hand, the compatibility relation can be written as 
(1 - a)m 2 



3R 2 



1 + R 2 + 1 + Rla 4 ^ 



with the normalized phase-speed, c, given as 



1 — a 



Rq 



cos 7. 



(lib) 

(11c) 

(lid) 

(lie) 



(12) 



(13) 



(14) 



From the next higher-order in e, we get the second-order approximations for variables which 
are of the following shapes 



(-4 +*)"?' +*(*+ *K >+ 



'(*+*K ) +'»(4+*)«'f , =o. 



(15a) 



(-i+*)^ (1) +K* + *)^ ) -^ (2) 



o, 



(15b) 



(-4 + A)«? ) +i(* + A)^+fi> u p ) = o, 



(15c) 



c (-* + *)^ (2)+m (* + *)^ ) = ' 



(15d) 



(2) 



3i? 2 



1 + i? 2 + a -V3^ 1 + i2 2 a 4/3 



(2) 



yielding the following second-order approximates 



(2) 



cfc 
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3i?o 2 [ v / TT^+a- 1 /Vl + ^« 4/3 ] M 2) (e,r) + ^ (2) (^r)] , 
[%>< 2 > (£, r) - d v ^ (v,r)]-l (£, r) + (rj, r)] , 

f ] = § [d^ i2) (M- V 2) (v,r)]+± [%y?« (e> r) + ^(D (r? , r) ] , 



(2) 

i 

,(2) 



(2) 



-3i? n 



(15e) 



(16) 



ra(l-a) 

where, v?^(£> r) and (p^(rj : r) describe the first-order amplitude evolution and v?^(£> T ) an d 
^ 2 \i], t) describe the second-order amplitude evolution components of two distinct solitary 
excitations in the oblique directions r/± and £j_ (r]± = — respectively In what follows, we 
will always employ the notations ip^ and ip^ instead of <£>^(£, t) an d V^HVi r ) f° r clarity. 



III. SOLITON DYNAMICS AND COLLISION PARAMETERS 



Proceeding to the next higher-order approximation, again by solving the coupled differ- 
ential equations in this approximation level and by making use of compatibility relation (Eq. 



8 



[T3T) and the previous plasma approximations, we obtain 



n. 



(3) _ N 



4 



Po(^T) + t/^ 



9y> 



5 



(i) 



S^ 3 



(1)1 



t _ N 
S 4 



(1) 



■ n^ , gave 



8 



1 



e 1 r.-W 



(i) 



a^ 3 



77+ 



Qo(£,r)-§ /^de 



a (!) 



9rj 



21 



4 



a 2 *f a 2 ^ 1 ' 



a? 2 



(17) 



+ 



where, the entities -F(£, r) and G(j], r) are assumed to be the homogenous solutions of 
differential equations. The coefficients in Eq. (|17p are as follows 

3/2 

A = u v t= . — — , ; — = cos 7, (18a) 



4^(1 - a)Ro 



+ a-V3^l + J R2 a 4/3 



F 3 , (u 2 + (1 - a)sin 2 7 ) 
6>/3(l - a)uj 2 



1 — a 



-i 3/2 



AT = 6v^(l - a)Rl 



1 -3/2 



COS 7, 



1 — a 



sec 7, 



(18b) 



(18c) 



E 1 = -3i? - 4 (1 + 2Rl{\ - a" 1 / 3 ) - a" 5 / 3 ) 



;isd) 



E> 



12R 2 



1 + R 2 + a- 1/3 yi + i2ga 4 /3 j , (18e) 

As it is realized already, the two first terms in Eqs. ( fl71) are secular and they diverge at 
£ ±oo and rj — > ±oo, and therefore must vanish. Thus, we obtain two distinct KdV 
evolution equations one for each traveling solitary structure. Furthermore, the next two 
terms in Eq. (TlTI) may become secular 29, [31] at the next higher-order and they must 



also vanish. The later conditions determine the collision parameters introduced in Eqs. 
(|E]). However, the full determination soliton dynamics is given by the following coupled 
differential equations 



dr 



+ M 



0. 



(19) 



Po(v,r) 



E 2 



(20) 



9 



(21) 



Qo&T) = ^Jtpfdt, (22) 



In order to get single-soliton solutions for Eqs. (1191) and (1211) . with multi-soliton solutions, 
we require that the perturbed potential components and their derivatives vanish at infinity, 
i.e. 



c hm o {^ 1) ,^,^-} = 0, C = t,V- (23) 



Then, we obtain 



ft - cosh2( i^oi)' 

tn — 3tt £Q A — ( ^B\\ 



(24) 



(25) 



,.(1) <Pr,0 

COsh ( ) 

^ = ^, A, = 

where, y?o and A represent the soliton amplitude and width, respectively, and uq is the 
relativistic Mach-speed. 

Close inspection of Eqs. ( fl~8l) . reveals that the KdV coefficients A and B in Eqs. ( fi~9l) 
and (I2T1) change sign at 7 = 7r/2, a critical value which determines weather the solitons are 
compressive or rarefactive. Moreover, it is remarked that the soliton amplitude is indepen- 
dent of the strength of magnetic field, while it strictly depends on applied field direction, 
cos 7. However, the soliton width depends on both magnitude and direction of the external 

nn 

magnetic field. This is in agreement with previous findings |1|. 1321] . 

The collision phase-shifts of solitary excitations are obtained using Eqs. ( 1201) and ( 122]) 
together with the KdV solutions (Eqs. ( 1241 and ( 1251) ) as 

P (r / ,r) = f^ A,tanh(^fl), (26) 
Qo(£,r) = §^oA ? tanh(^«. (27) 
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We can write, up to order 0(e 2 ) 

(28) 



£ = + / y + m2 + ct) -e 2 g^ A r ,tanh(^^) + 0(e 3 ), 



77 = + /«/ + mz- ct)] - e 2 |^^ A ? tanh(^-^) + 0{e 3 ), 

We can also calculate the overall phase-shifts by comparing the phases of each wave long 
before and after the collision event in the following form 

a p p _ p 

1 post— collision 1 past— collision 

lim [e(kx + ly + mz + ct)] — lim [e(kx + ly + mz + ct)], 

£,=0,n— >+oo £=0,n— >-oo (29) 
AQq Q post— collision Q past— collision 

lim [e(kx + ly + mz — ct)] — lim [e{kx + ly + mz — ct)], 

The phase parameters, APo and AQo; therefore, present the overall phase-shifts of solitary 
structures labeled "51" and "52" in the head-on collision. Finally, making use of Eqs. ( |26f) . 
fl27|) and (jSJ), we obtain the following simplified expressions 



AP = £ 



l+2ij f 2 ,(l-a- 1 / 3 )-a- 5 / 3 



/3 



AQo = -e 



2 J R2( A /T+fi|+a-l/3 % /i +jR 2 Q 4 
l+2 J Rg(l-q- 1 / 3 )-a- 5 / 3 



2R2( A /I + -^ +Q -l/3 A y 1+iJ 2 Q 4/3) 



^rjoA^, 

(30) 



IV. NUMERICAL ANALYSIS 

In this section we evaluate the nonlinear wave dynamics using typical values of mass- 
density for a degenerate white dwarf. Figure 1 depicts the variation of amplitude and width 
of a solitary excitations with respect to various parameters in a magnetized degenerate 
plasma with relativistically degenerate electrons and positrons. It is remarked from Figs. 
1(a) and 1(b) that, the increase of the mass-density enhances both the amplitude and width 
of the excitations for fixed other plasma parameters. It is also indicated that the only 
compressive solitons are present in this plasma for all given parameters. However, regarding 
the relative positron density variations, one observes that, the soliton amplitude/width 
increases/decreases as the relative positron density increases. Moreover, as indicated by 
Figs 1(c), the amplitude of soliton is minimum when the magnetic field is parallel to the 
propagation direction and increases as the the field direction approaches to the normal, 
7 = 7r/2 position. The soliton width, on the other hand, shows (Fig. 1(d)) a maximum 
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value at field-angles denned by 



1 



2 1 a 



arccos 




1 — a 



(31) 



where the maximum value of this quantity itself is ^ wm ~ 54.73°. This is in complete 
agreement with the cases presented in Ref. |]J treated with a classical model for ultra- 
relativistic and ordinary electron-positron-ion plasmas. It is from Eq. (I3"T|) remarked that 
this maximum width field-angle is not affected by the mass-density of plasma. However, 
the value of the maximum width increases as the plasma gets denser. Figures 1(e) and 
1(f) reveal that, while the amplitude increases monotonically with increases in the relative 
positron number-density, the width has a maximum value which enhances and tend to lower 
a- values as the plasma mass-density increases. 

Figure 2 shows variations of the head-on collision phase-shift for interacting solitons. 
It is observed that, the head-on collision phase-shift is always positive indicating that the 
post-collision parts of solitons always move ahead of initial trajectories [33j. The value of 
collision phase-shift increases with increase of the plasma mass-density, as it is concluded 
from Fig. 2(a). It is also observed from Fig. 2(b) that this parameter is increased as the 
magnetic field direction moves from the direction parallel to propagation to normal position, 
however, this effect gets more pronounced when the mass-density of plasma is higher. On 
the other hand, the variation of collisional phase-shift with respect to the relative positron 
concentration, presented in Fig. 2(c) is much similar to that of the width of soliton shown 
in Fig. 1(f). Furthermore, Fig. 2(d) indicates that the increase in the relative magnitude of 
magnetic field, u, has also significant effect on the collision of solitons, so that, the higher 
is the field strength the lower is the shift in the soliton trajectory in the head-on collisions. 

V. CONCLUDING REMARKS 

Using the extended Poincare-Lighthill-Kuo (PLK) reductive perturbation method we 
studied the oblique propagation and quasi-elastic head-on collisions of solitary structures 
in a fully degenerate magnetized electron-positron-ion plasma with relativistic ultra-cold 
electrons and positrons in the framework of quantum hydrodynamics model. It is observed 
that, plasma parameters such as mass-density, normalized magnetic field strength, its an- 
gle with respect to the soliton propagation and the relative positron number-density play 
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important roles on the propagation as well as collision dynamics of solitary structures in a 
Fermi-Dirac plasma. Current findings can elucidate new aspects of nonlinear dynamics in 
dense astrophysical objects such as white dwarfs and may as well be applicable to 
the laboratory produced inert ial- confined fusion plasmas. 

Appendix A: Normalized equations in strained coordinate 



r 2 drij _ „dni _ r 2„ dQo drij , p §n± , r 2 r dPp drij , j drijUj , 
c dr u d£ c u 9? dr) &n 9r) d£ ^ d£ 

r 2i„ dQo dnjUj _|_ i drijUj , 2i„ dPo drijUj , i drijVj , ^.2j dQo dnjVj , 

fc ^ a? a?? ^ ^ a?? ^ fc ^ a?? ai ^ 1 a? ^ fc 1 a£ ar? " r 
a»? dr/ d£, a§ a£ a^ ar; 

c 2 dPodriiUH , _ n 

L " L dn di i- • • • — u, 

r 2 a«, _ a^ _ 2 ago a«, , au, , 2 a_Po a^j , k a«j . 
fc ar c ae fc c a? a»? ^ a»? fc c a»? a? ^ a? ~ r 

t r*Ui d £ dv -1- i«i ^ -1- t 9j? ^ -1- if i ^ -1- t tfj ^ 9?) -i- 

Z^ff + £ 2 ^f^ + rm^ + £ 2 mw .^o ^ + mWl ^+ 

^2 ap a M j 1 r,a£ , 2j,aQo a^ , > a^ , 2udPo 9jp _ _i_ - n 

r 2 _ ^a^ _ -.2 ago 1 a^ , 2 ap a^ , /_ a^ . 
c ar a? c a? a?j T ar? T c a?? a? T ac T 

r 2u dQodvi 1 U-.dvi _|_ F 2h, 1 .dPodvi , /^.a^ , P 2 ] „ , dQo 9ih _i_ 

t AiUj a ^ 9?? -r ^ -r t ^ a? -r tf ? a ^ -r t a ^ a?? -r 

t + ^ + + 1% + + ^ + • • • = 0, 

r 2 dwj _ dwi_ _ p 2 dQo dwj , dm , 2 dPo&m i I... a^j , 

fc ar u a? c ^ a? a^ ^ ^ a?? ^ c ^ a?? ac ^" a^ T 

r 2; dQ dwj , U^.dm i p 2l. 1 .dPodwi , U.dvH , P 2u.dQo_dwi\ 

fc a? a?; ^ a?; ^ fc ar, a? ^ ty * a? ^ fc tt7 « ac a^ ^ 



e 2 



"av 








dr) 2 


-{ 



^ [(1 + 2Rl) - a-^{\ + 2Rla^)] } 



+ ... = 0. 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



13 



[1] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 092304(2010). 

[2] F. Haas, L. G. Garcia, J. Goedert, and G. Manfredi, Phys. Plasmas 10, 3858(2003). 
[3] Shukla P K, Eliasson B, "Nonlinear aspects of quantum plasma physics" Phys. Usp. 51 
53(2010). 

[4] G. Manfredi, Fields Inst. Commun. 46 263(2005). 

[5] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 052302(2010). 

[6] M. Bonitz, D. Semkat, A. Filinov, V. Golubnychyi, D. Kremp, D.O. Gericke, M.S. Murillo, 
V. Filinov, V. Fortov, W. Hoyer, J. Phys. A, 36 5921(2003). 

[7] S. Chandrasekhare, Science, 226 4674(1984). 

[8] A. A. Mamun and P. K. Shukla, Phys. Lett. A, 374, 4238(2010). 

[9] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 072101(2010). 
[10] M. Akbari-Moghanjoughi, Astrophys. Space Sci., (2010) DOI 10.1007/sl0509-010-0499-3. 
[11] N. J. Zabusky and M. D. Kruskal Phys. Rev. Lett. 15 240(1965). 
[12] H. Washimi and T. Taniuti, Phys. Rev. Lett. 17 996(1966). 
[13] M. Oikawa and N. Yajima, J. Phys. Rev. Sco. Jpn. 37 1093(1973). 

[14] Shapiro, S. L. and Teukolsky, S. A., Black Holes, White Dwarfs and Neutron Stars, (Wiley, 
New York), (1983). 

[15] M. J. Rees, in: G. B. Gibbons, S. W. Hawking, S. Siklas (Eds.), The Very Early Universe, 

Cambridge Univ. Press, Cambridge, (1983). 
[16] H. R. Miller, P. J. Witta, Active Galetic Nuclei, Springer- Verlag, Berlin, p. 202, (1987). 
[17] P. Goldreich, W. H. Julian, Astrophys. J., 157, 869(1969). 
[18] F. C. Michel, Rev. Mod. Phys. 54, 1(1982). 

[19] L. O. Silva, R. Bingham, JM. Dawson, J. T.Mendonca, and P. K. Shukla, Phys. Rev. Lett. 
83 2703(1999). 

[20] C. M. Surko, M. Levelhal, W. S. Crane, A. Passne, F. Wysocki, T. J. Murphy, J. Strachan, 

and W.L. Rowan, Rev. Sci. Instrum. 57 1862(1986). 
[21] C. M. Surko, T. Murphy, Phys. Fluids, B2 1372(1990). 
[22] R. G. Greeves and C. M. Surko, Phys. Rev. Lett. 75 3847(1995). 
[23] V. I. Berezhiani, D. D. Tskhakaya and P. K. Shukla, Phys. Rev. A, 46 6608(1992). 



14 



[24] S. Chandrasekhar, "An Introduction to the Study of Stellar Structure", Chicago, 111. , (The 

University of Chicago press), (1939) p. 360. 
[25] A. Rasheed, G. Murtaza and N. L. Tsintsadze, Phys. Rev. E., 82 016403(2010). 
[26] Nam C. Lee, Phys. Plasmas, 16 042316(2009). 

[27] Mouloud Tribeche and Soufiane Boukhalfa, Astrophys. Space Sci., (2010) DOI: 

10.1007/sl0509-010-0549-x 
[28] A. Jeffery, T. Kawahawa, Asymptotic Method in Monlinear Wave Theory, (Pitman, London, 

1982). 

[29] Masayuki OIKAWA and Nobuo YAJIMA, Phys. Soc. Jpn., 34, 1093(1973). 
[30] E. Infeld, G. Rowlands, Nonlinear Waves, Solitons and Chaos, Cambridge Univ. Press, Cam- 
bridge, (1990). 

[31] Nobumasa Sugimoto and Tsunehiko Kakutani, J. Phys. Soc. Jpn. 43, 1469(1977). 
[32] A. Esfandyari-Kalejahi, M. Mehdipoor and M. Akbari-Moghanjoughi, Phys. Plasmas, 16 
052309(2009). 

[33] M. Akbari-Moghanjoughi, Phys. Lett. A, 374, 1721(2010). 



15 



FIGURE CAPTIONS 



Figure 1 

(Color online) Variations of soliton amplitude (left column of figure) and width (right 
column of figure) with respect to magnetic field direction 7, relative positron concentrations, 
a, and normalized mass-density of plasma. The dash sizes in all plots are appropriately 
related to the values of the varied parameter. The values of u^ = u V; q = 0.1 and e = 0.1 
are used for all plots. 

Figure 2 

(Color online) Variations of the head-on collision phase-shift with respect to with respect 
to magnetic field direction 7, relative positron concentrations, a, and normalized mass- 
density of plasma, when other parameters are fixed. The values of u^q = u r),o = 0.1 and 
e = 0.1 are used for all plots. The dash sizes in all plots are appropriately related to the 
values of the varied parameter. 
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